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ABSTRACT 

A Nusselt-type analysis was performed for laminar film 
condensation on the inside of a rotating truncated cone 
with small half cone angles. This analysis included the 
interfacial shear between the vapor and condensate, the 
vapor pressure drop, the thermal resistance in the condenser 
wall, and the condenser outside cooling mechanism. An 
approximation of the analytical model made it possible to 
find a numerical solution for small half cone angles greater 
than zero. 

A non-capillary rotating heat pipe containing an 
evaporator, condenser, and distilled water as the working 
fluid was tested. It was rotated at 702 and 1404 RPM, and 
the heat transfer rates of the heat pipe were determined 
experimentally for different saturation temperatures cor- 
responding to electrical power inputs ranging from 1 kW 
to 9 kW. 

The experimental results showed that the non-capillary 
rotating heat pipe was an effective heat transfer device. 

The approximate numerical solution conservatively predicted 
the heat transfer rate with a deviation of 18% at 702 RPM 
and 5% at 1404 RPM. 
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I. 



INTRODUCTION 



A. CONVENTIONAL HEAT PIPF 

A heat pipe is a closed container capable of transfer- 
ring large amounts of heat at nearly isothermal conditions. 
Conventionally it consists of four main parts: (1) evap- 

orator, (2) condenser, (3) working fluid, and (4) capillary 
wicking structure to force the condensate back to the 
evaporator. The working fluid is evaporated in the evap- 
orator section of the heat pipe by the heat source. As the 
vapor pressure increases, the vapor flows toward the con- 
denser end, transporting the heat as the latent heat of 
vaporization. At the condenser end, the heat is removed 
from the vapor as it condenses. The condensate is then 
returned to the evaporator by the capillary pumping action 
within the wicking structure. 

Studies of the heat transfer capabilities have been made 
on the conventional heat pipe [6, 10]* with the conclusion 
that, for ordinary fluids and moderate heat fluxes, the heat 
pipe is limited in its capability to transfer heat by the 
rate of condensate that can be returned to the evaporator. 

B. THE ROTATING, NON-CAPILLARY HEAT PIPF 

The rate of condensate flow may be increased by removing 
the wicking structure from the condenser. If the condenser 

*Numbers in brackets indicate references listed in the 
bibliography. 
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end is then internally tapered and the heat pipe is rotated 
about its longitudinal axis, a centrifugal force field is 
set up which will overcome both the viscous losses in the 
system and the gravitational field. 

Ballback's study [10] of the rotating heat pipe indi- 
cated that its heat transfer capability is limited by the 
amount of heat which can be removed from the condenser. In 
his study, Ballback performed a Nusselt-type analysis for 
the condensation on the inside of a truncated cone. In 
order to obtain a closed form approximate solution for the 

film thickness, 6, he made the following assumptions: (1) 

d 6 

« tan <J> (large cone angles) , (2) no shear stresses between 

the vapor and condensate, and (3) no thermal resistance in 
the condenser wall. Daley [11] was able to improve on the 
theory introduced by Ballback by taking into consideration 
the thermal resistance in the condenser wall and the outside 
cooling mechanism. Furthermore, Daley neglected the assump- 
tion that ~ « tan <f>. Using the velocity profile for the 
condensate developed by Ballback, he arrived at a non-linear 
second order differential equation for the film thickness , 

6 . Daley numerically integrated this expression using a 
Runge-Kutta-Gill numerical integration scheme with an 
IBM 360 MOD 67 digital computer. To start the integration, 
he used the conditions that <5 = 6i and = tan<J> at X = 0 . 

It was postulated that 6i was a function of the over-fall 
condition at the step located at the condenser exit. Daley 
arrived at an expression for the minimum film thickness, 6m, 
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at the condenser exit, using the approximation that the flow 
over the step approximated the free fall condition for open 
channel flow. Having found the film thickness from the 
second, order differential equation at X = Lc, the condenser 
exit, Daley compared this value with 6m at Lc. When the 
difference was less than .004, a solution was found. Using 
this technique, Daley was able to find the heat flow out of 
the condenser as a function of RPM for 0° half cone angle 
only. As the cone angle was increased, however, instability 
of the highly non-linear second order differential equation 
for the condensate film thickness caused the equation to be 
extremely sensitive to the initial value of the film thick- 
ness . Consequently, no solutions were found for the conden- 
sate film thickness for half cone angles greater than zero. 

C. THESIS OBJECTIVES 

The objectives of this study were to: (1) continue the 

theoretical study of the heat transfer capability of the 
rotating heat pipe and obtain solutions for <f> > 0° , (2) 

obtain experimental data from the heat pipe designed and 
assembled during Daley's study, and (3) compare the analyti- 
cal and experimental results. 
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II. THEORETICAL PROGRAM 



In order to improve on the analysis performed by 
Ballback and Daley the effects of vapor pressure drop and 
interfacial shear stress between the liquid and the vapor 
were included. Furthermore, the thermal resistance in the 
condenser wall and the outside cooling mechanism were con- 
sidered. Using the coordinate system in Figure 1, for 
small half cone angles, momentum, continuity, and energy 
equations were developed for an infinitesimal fluid element. 
In addition, assuming one dimensional turbulent flow in the 
vapor, momentum and continuity equations were developed for 
the isothermal vapor flow. 




Figure 1. Coordinate system and geometry 
for condensate film analysis. 
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A. CONDENSATE MOMENTUM EQUATION (X~Direction ) 

By assuming that the momentum changes in the condensate 
are small a static force balance may be taken on an infini- 
tesimal liquid element in the x-direction as shown in 
Figure 2 . 

0 = " §§ + p f“ 2 r sin * (1 > 

where 

T = shear stress 
p = pressure 

x = co-ordinate measuring distance along surface 

y = co-ordinate measuring distance normal to surface 

r = radius 

<f> = half cone angle 

p£ = density of liquid 

0 ) = angular velocity 

B. CONDENSATE MOMENTUM EQUATION ( Y-Direction) 

Similarly, using Figure 2, a force balance in the 

Y-direction may be taken 

0 = - r cos-<J> (2) 

Making use of Figure 2, the radius at any point in the fluid 
condensate may be expressed as: 

r(x,y) = R q + x sin <f> - y cos <f> (3) 

where 

R = minimum wall radius in the condenser section 
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Figure 2. Cross-section of infinitesimal liquid 
element . 
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Substituting equation (3) for r(x,y) into (2), and inte- 
grating (2) from p to p v and y to 6 yields: 

2 

p(x,y) = P v (x) + p f u) cos <|>{(R + x sin<f>) (<$-y) 

- j cos <K6 2 -y 2 ) } (4 ) 

Differentiating equation (4) with respect to x gives 

3p dp v 2 , //T _ , « . p. dfi 

9 ^ = + p f w cos 4){(R o + x sin cf>- 6 cos <J>) ^ 

+ (6-y) sin c{>} (5) 

Substituting equation (5) into equation (1) and integrating 
from y to 6 and t to -t gives 

dp 7 

t = - (6-y) - T v “ P f w cos <f>{(R Q + x sinc{)-6cos (f>) (6-y)} 

2 

+ p^co sin <J> {R q +xsin <J> - 6 cos <j> }(6-y ) 



9u (x,y) 



f 3y 



(6) 



where 

y f = viscosity of the liquid 

u(x,y) = velocity of the liquid 

Integrating equation (6) from 0 to y and from 0 to u gives 
an expression for the velocity of the fluid: 



, , i dp v ,, zi T v y 

u(x,y) = - — 5T" (6y “ ” — 



P f W* 



J f L 



R q + x sine} 



• A A d6 

sm(J)-cos(j)^- 



2 ”) 

6v - \\ 



(7) 
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This expression for u(x,y) agrees with Ballback's except 
for the shear stress term. 



C. CONTINUITY EQUATION 

It is apparent from Figure 3 that due to the conserva- 
tion of mass, with condensation 



or 



dw dw 
f _ v 

dx dx 



w 

v 



( 10 ) 



(10a) 



where 

Wf = mass flow rate of fluid 
w^ = mass flow rate of vapor 

The mass flow rate of the liquid may be written as 
6 

w f = f P f u 2 7T r dy (11) 

o 

The integrating of equation (11) using equations (3) and 
(7) yields 



'f = 2 ^pf { - ^ 



dp 

dx 



[(R + x sin 0) 



_5_ 

24 



6 4 ] } 



T r 2 p 3 

+ 27rp f {- —■[ (R q + X sin <f>) “2 cos <f> } 



2 2 

27rp f a) 



d6 x 



(R q + x sin cf) - <5 coscj))(sin c{) - cos ^-) 



[ (R q + x sin - cos * 6 4 ] 



(12) 



The mass i ate of flow of the vapor may also be written 
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Figure 3 . 



Cross section of ^ 
liquid and vapor elements. 
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( 13 ) 




where 



R = R + x sin <b - 6cos <f> 
o r Y 



D. ENERGY EQUATION IN CONDENSATE 

From Figure 2, using an energy balance, the rate of 
condensation may be written 



T g = saturation temperature of vapor 

T m = inlet coolant temperature 

t = thickness of condenser wall 

6 = film thickness of condensate 

k^ = thermal conductivity of condensate 

k^ = thermal conductivity of condenser wall 

h = outside condenser wall heat transfer 



coefficient , 

E. VAPOR MOMENTUM EQUATION 

Using Figure 3, and neglecting the momentum of the 
condensing vapor, a force balance may be taken in the 
x-direction. 




2ir (R q + x sin (j> ) (Tg-TJ 



(14) 



dx fg dx 




where 



P V A " (p v +dp v ) (A-dA) - (p v +-jX)dA - T v dA w 



(w v +dw v ) (v+dv) - w v v 



(15) 



where 
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P v -- pressure of vapor 
v = velocity of vapor 

A = cross sectional area of vapor element 

A w = surface area over which the shear 
stress acts 

Neglecting higher order terms, and using the fact that 
dA^ = 2‘irRdx equation (15) reduces to 



dx 



P v 
v 



dv 

dx 




_d_ 

dx 



(w v ) + 2 



R 



(16) 



1 2 

Using the expression T v = - p v v f, where f = friction factor, 
and substituting into equation (16) yields 



dp 

dx 



dv 

dx 






+ 



p v 2 f 
K v 

R 



(17) 



Expanding the derivative in equation (17) yields 



dp 

dx 



2P v v 



dv 

dx 



2v 2 p 



v dR 
dx 



p v 2 f 
v 

R 



(17a) 



The first two terms in equation (17) represent the momentum 
effects of the vapor while' the last term represents the 
usual friction effect. For the low condensate velocities 
which were expected during the experimental runs , according 
to Bergelin et al. [2] the friction factor is approximately 
equal to the friction factor used in smooth tubes (that is , 
the effect of the liquid film on the condenser surface is 
negligible) . According to Rohsenow and Choi [3] , the 
friction factor for turbulent flow in a smooth tube for vapor 
Reynolds numbers between 2000 and 50,000 is expressed as: 
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f = 



.0791^ 



(Re) 35 = 




(18) 



F. DETERMINATION OF FILM THICKNESS AND HEAT TRANSFER RATE 

Equations (10) , (12) , (13) , (14) , and (17a) can be 

further reduced to two simultaneous ordinary first order 

differential equations. Using equation (14) with the use 

dwf 

of equations (10) and (13) for ^—yields: 



dv 

dx 



2 (R q + x sin 4>) (T 




T ) 

CO' 




2v 

R 



sin 4 > - cos 4 > 



d_6_*1 

dxj 



(19) 



At the same time, using equation (12) in conjunction with 
(10a) and (13) we obtain 



P v ttR“v = 2Trp f j 



f 1 dp v . , x 6 3 *5*4, 

— UK + x cos ^ 24 u J 



+ 2ttp, 



y f dx 

T r 2 

-r~ [ (R + x sin4>) 

M jZ ^ 



2 " - cos 4> 3 -] 



27Tp f 2 0) 2 d6 

+ (R + x sin (J) - 6cos4>) (sin 4)- cos — ) 






. 6 . 5 r 4 

. [ (R q + x sm4>) — - cos 4) 24 - 6 



( 20 ) 



If equation (17a) is now substituted into equation (20) , 



equation for is obtained: 



d6 

— = tan c 
dx 



2 


[Pi p 2 ' 




-V f R 


0n dv 
2p l dx 


pvt 


R + 2 


+ P v v 


2 P £ 


P 1 cos 4^ 


2 

Pf to R + ■ 


2p v v 2- 

R 





( 21 ) 



where 



5 * 4 
24 6 



6 

= (R 0 + xsin4>) ~ 2 ~ - cos 

<$ 2 6 3 

p 2 = (R q + xsin4>) cos 4> — 
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Equations (19) and (21) are two simultaneous first 
order ordinary differential equations with unknowns 6 (x) 
and v(x). A Runge-Kutta numerical integration was used to 
solve them with the IBM digital computer. It was hoped 
that by using first order differential equations rather 
than a second order, as Daley did, that the sensitivity to 
initial 6 would be reduced. Figure 4 is a flow chart of 
the program used to solve the differential equations. 

The integration was started by assuming initial con- 
ditions, 6=6^ and v=0 at x=0. With these initial conditions 
= tan<f) at x=0. Figure 5 shows a schematic representing 
the calculated 6 (x) for different initial values of <$ . 

Note that if 6. is too small the calculated 6 becomes 
negative before the condenser overfall (physically impossible 
solution) . If 6^ is too large the calculated 6 grows at a 
constant slope equal to tan<f> (physically unrealistic). With 
the proper choice of 6^ the calculated 6 reaches a maximum 
along the length of the condenser and then tapers off to a 
value close to zero at the overfall. The exact value of 6 
at the overfall is not known. It was observed however that 
the calculated value of 6 at the overfall was approximately 

.1 of <$ . , and the heat transfer rate was insensitive to an 
i 

exact value at the overfall. This result agrees with Nimmo 
and Leppert [8] who experimentally found for laminar film 
condensation on a finite’ horizontal surface that the heat 
transfer rates were practically independent of the end 
conditions . 
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Figure 4. Flow chart of computer program. 
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tan <p - 



p v v 2f 




+ p V 
v 


r p f R dvi 

2p f 2p l dx J 


P 1 COS (f> 


0 2p v^ 

[pf“ 2 » + -i- 





dv 

dx 



2 (R q + x sin <f>) (T 



TJ 



r 6 t 11 


2 


p vk7 + r + r 

|_ f w 


R h fg 



2 v sincp 
R 



( 22 ) 



(23) 



Equation (22) may be further simplified by an order of mag- 
nitude analysis. Neglecting the second terms of and P 2 

and using the expression, R = R + x sin6 with the 

o 

previously observed result that p n ^-<< -- — 

1 dx 2p, 



dv „ p f R2 „ Pi P 2 



and — « — , 

equation (22) may be rearranged in terms of a cubic equation 
for 6 (x) : 



R sind 



2p v V 



Pf“ R + R 

P f R 2 
■ p v v 2 ^~ = 



p v 2 f 
K v 






(24) 



Equations (23) and (24) were solved simultaneously for 
v(x) and 6 (x) . Equation (-23) was numerically integrated 
using 1000 steps of integration over the cooling length of 
the condenser. The integration was started with initial 
conditions 6^ = 0 and v = 0 at x = 0. Having found v(x) at a 
particular step, equation (24) was solved by a subroutine. 
This value of 6 was used to find v(x) for the next step of 
integration. The results of approximation 1 are recorded 
in Table I for half cone angles of .1°, .2°, and .3°. 

Using approximation 2 on equation (21) , equation (24) 
was further simplified and rearranged to yield the following 
expression for 6 (x) : 
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TABLE I 



THE HEAT FLOW, Q (Btu/hr) , FOR 
DIFFERENT HALF CONE ANGLES 



<f> 


NUMERICAL 
SOLUTION 
Eqs . (19) & (21) 


APPROXIMATE 
SOLUTION 1 
Eqs. (23) & (24) 


APPROXIMATE 
SOLUTION 2 
Eqs. (23) & (25) 


0° 


8,505 


- 


- 


.1° 


10,082 


9,836 


10,290 


.2° 


11,173 


11,033 


11,332 


.3° 


11,904* 


11,769 


11,998 



6(x) 



3p p .v 
f 

2 2 . 
2p f a) sm 




(25) 



Using equations (23) and (25) a numerical solution was found 
for v(x) and 6 (x) in the same manner as the solution for 
approximation 1. The results are tabulated in Table I for 
half cone angles of .1°, .2°, and .3°. 

Table I shows that approximation 1 is lower than 2. 

This is to be expected since approximation 1 includes both 
interfacial shear between the' vapor and condensate and the 
pressure drop in the vapor. At the same time approximation 
1 is less than the more accurate solution of equations (19) 
and (21) , thus giving a conservative result with a deviation 



*In finding the solution for .3°, 10,000 steps were 
used in the integration. This was necessary because as the 
half cone angle became larger than .2°, the numerical inte- 
gration became very sensitive to the step size. 
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of only 2.5%. Consequently, approximation 1 v/as used to 
arrive at a theoretical solution for 6 (x) and Q for the 
one degree half cone angle condenser used in the experi- 
mental program. The results are plotted in Figures 9 
and 10. 
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II. EXPERIMENTAL PROGRAM 



A. DESCRIPTION OF EQUIPMENT 

In order to facilitate the description of the entire 
experimental system, the system is sub-grouped into four 
parts: (1) evaporator, (2) condenser, (3) instrumentation, 

and (4) auxiliary equipment. Included in the description 
are drawings which are intended to aid the reader in 
visualizing the different components of the system. 

Figure 6 is a schematic representation of the entire system 
including the auxiliary equipment. Figure 7 is a cross- 
sectional drawing of the assembled heat pipe. Figure 8 
is a photograph of the heat pipe and associated equipment. 

1 . Evaporator 

The evaporator, modeled to approximate the rotating 
boiler tested at Lewis Research Center [17] , is made of 304 
stainless steel. It has an inside diameter of 3.125 inches, 
an outside diameter of 3.6*35 inches, an exit diameter of 
2.625 inches, and an overall length of 5.905 inches. 

The actual heated length is 3-7/8 inches around which 
an 11-gauge Chromel-A heater wire in a 3/16 inch inconel 
sheath is helically wound. The procedure of bonding the 
heater wire to the evaporator is an exacting process in 
order to ensure a uniform heat flux axially along the heater 
surface and to guard against electrical shorting between 
the heater wire and other components. The procedure was as 
follows : 
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FIGURE SCHEMATIC DIAGRAM OF APPARATUS 
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FIGURE 7 CROSS SECTION OF ROTATING HEAT PIPE 
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Figure 8. Photograph of heat pipe and associated equipment. 




a. A sheet of gold-nickel foil, .003 inch thick, 
was tack welded to the outside evaporator surface. 

b. The heater wire was helically wound over the 
gold-nickel sheet with two wraps of .020 inch diameter gold- 
nickel wire wound intermittently between the heater coils. 

c. The assembly was then brazed in a hydrogen 
atmospherically-controlled oven . 

d. Sauereisen cement was baked around the outside 
of the heater wire using the heater wire as the heat source. 

e. Asbestos insulation was wrapped around the 
heater wire and tied off with glass tape. 

f. Fish paper was placed into the heat sink on the 
condenser end of the evaporator to insulate the heater 
coils from the evaporator. 

Power was furnished to the heater coils by a DC motor- 
generator capable of delivering 150 amperes at 220 volts. 

The electrical power was passed by a brush assembly con- 
sisting of eight spring-loaded graphite brushes through two 
bronze collector rings to the helically wound heater wire. 
The rings were supported by twelve 1/4 inch stud bolts. 
Phenolic spacers were used to align the collector rings and 
to insulate the bolts from the rings. 

In order to reduce the axial heat losses, two radial 
1/16 inch wide grooves were machined in the outside surface 
to within 1/16 inch of the inside boiler wall at either end 
of the heating element. 



30 



The outside face of the evaporator is closed by the 
insertion of two pyrex glass view ports, 3-1/2 inches in 
diameter. The inner view port rests on a neoprene O-ring 
seal and is separated from the outer view port by a 
3/32 inch compressed fiber gasket. The view ports are held 
in place by a stainless steel end flange secured in place 
by twelve 1/4 inch stud bolts. The original design of two 
solid glass view ports was changed by fusing the inner view 
port to a glass tube which passes through the outer view 
port. This tube allows a vacuum to be drawn from the inside 
of the heat pipe during the filling procedure, and is easily 
sealed . 

The inside diameter of the evaporator is stepped 1/2 
inch from the exit diameter to allow for the formation of a 
liquid annulus in the evaporator. 

2 . Condenser 

The condenser section is machined of stainless steel 
in the shape of a truncated cone. The half cone angle is 
1°. The base of the condenser mates to the evaporator exit 
and has an inside diameter of 1.809 inches and tapers over 
a 10 inch length to a diameter of 1.460 inches at the 
truncated end. At this point, the condenser extends an 
additional 1-1/2 inches to house a cylindrical end plug. 

A step is machined at the evaporator-condenser inter- 
section where the condensate flows from the condenser to the 
recessed inside diameter of the evaporator. The first inch 
of the condenser surface near the evaporator is assumed to 
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be adiabatic and 9 inches of the condenser's overall length 
is considered as the cooling surface. The wall thickness 
of the condenser is 1/16 inch in order to reduce its thermal 
resistance . 

The condenser section has two flanges. The large flange 
mates to the bearing seat flange on the evaporator and is 
held in place by the 1/4 inch stud bolts. The other flange 
mates to a flange on the end plug held together by eight 
1/8 inch stud bolts. The end plug has an extension on its 
open end on which are mounted four phenolic junction boards. 

3 . Instrumentation 

A schematic drawing of the instrumentation set-up 
is shown in Figure 6. 

a. Temperature Measurements. 

Temperature measurements are obtained from nine 
grounded copper-constantan thermocouples and two quartz 
thermometers. The thermometers are used to measure a 
temperature profile in the evaporator wall, the liquid 
annulus temperature in the evaporator, the vapor saturation 
temperature in the evaporator, the outside wall temperature 
of the condenser, and the temperature at the location of the 
pressure transducer. 

In order to measure the temperature gradient in the 
evaporator wall, four thermocouples with a 1/16 inch dia- 
meter stainless steel sheath are utilized. Each thermo- 
couple is placed within a 1/16 inch well drilled 3-1/2 
inches deep into the evaporator wall. The thermocouple 
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junction is coated with epoxy to give it enough strength so 
that it will not bend when placed against the end of the 
well. The epoxy is sanded such that the junction makes 
contact with the wall. 

One thermocouple is located 1/16 inch radially from the 
inside of the evaporator, another is located 1/16 inch 
radially in from the outside evaporator surface, and the 
remaining two are located opposite one another in the middle 
of the 1/4 inch wall at the same radial position to check 
on the uniformity of the temperature field. 

The leads from all of the thermocouples in the evapor- 
ator wall are brought through slots in the evaporator and 
forward condenser flanges, down the side of the condenser, 
and through the rear condenser flange to phenolic junction 
boards. From the boards, the leads go inside the shaft and 
down to the slip ring unit which is attached to the main 
assembly by a coupling mounted at the end of the shaft. 

A thermocouple is suspended in the liquid annulus to 
measure the liquid bath temperature. The remaining thermo- 
couple in the evaporator is suspended in the vapor space to 
measure the saturation temperature. 

Since these two thermocouples protrude into the inside 
of the heat pipe where a vacuum is held, precautions were 
taken to ensure that a seal was maintained. The ends of 
the two thermocouples were enclosed in stainless steel 
flanges and the leads were soldered at the slots on the 
condenser forward flange to keep air from passing through 
the slots into the evaporator. 
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In order to measure the outside condenser wall temper- 
ature, two thermocouples were soldered to grooves on the 
outside of the condenser at locations of 1/3 and 2/3 of its 
length respectively. 

The remaining thermocouple was used to measure the 
temperature at the pressure transducer. The junction is 
soldered to the inside of the pressure transducer arm 1/32 
inch from the surface. 

The temperature of the coolant is measured before and 
after spraying by quartz thermometers mounted in mixing 
tubes located above and below the condenser spray cooling 
cylinder . 

b. Pressure Measurements 

The saturation pressure of the vapor at the 
small end of the condenser is measured by a semi-conductor 
pressure transducer mounted in the vapor space. The trans- 
ducer is temperature compensated and is threaded on a 1/4 
inch diameter arm that extends from the shaft through the 
end plug and about three ' inches into the vapor space. The 
leads of the probe pass through the transducer arm up 
through the shaft to the junction boards. From there, the 
leads pass back through the shaft to the slip ring unit. 

c. Coolant Flow Rate Measurements 

The coolant flow rate is measured by passing 
the coolant through a flow rotameter before the coolant 
enters the spray nozzles. 
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d. Rotational Speed Measurements 

The rotational speed of the heat pipe is 
measured by two meters. A Beckman Eput Meter measures the 
RPM using a magnetic pickup, and an electronic strobe scope 
is used as a check. 

e. Electrical Power Measurements 

Power to the heater coils is measured by an 
ammeter and voltmeter connected in series and parallel 
respectively with the power lines returning from the heater. 
The power level is controlled by a field-rheostat placed in 
parallel with the DC motor-generator. A 150 ampere circuit 
breaker is included in the circuit. 

f. Slip Ring Unit 

The slip ring unit allows the rotating thermo- 
couple leads and pressure transducer leads to pass electric 
current to stationary leads. The slip ring unit used in 
the experiment utilizes liquid mercury as a viscous contact 
medium between the rotating and stationary leads. The 
frictional heating noise is less than 10 microvolts at 
speeds up to 4000 RPM. The slip ring unit has 22 terminals 
to accommodate the nine thermocouples and pressure trans- 
ducer. The output from the slip rings is fed to a Hewlett- 
Packard 2010C Data Acquisition System containing an inte- 
grating digital voltmeter, guarded data amplifier, and 
guarded cross-bar scanner. This system provides printed or 
visual output with an accuracy of ± .5 microvolts. 
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4 . Auxiliary Equipment 



The auxiliary equipment includes: (a) drive 

assembly, (b) test stand, (c) safety shields, and (d) spray 
assembly . 

a. Drive Assembly 

The drive assembly consists of a drive motor, 
end plug, shaft, pulley and bearing. The drive motor drives 
the shaft by means of a 2.65 inch diameter pulley and 
V-belt. The driving motor is a 2 horsepower, 3-phase 
variable speed motor capable of a maximum speed of 4000 RPM „ 
The speed is controlled by a pushbutton electric remote 
control unit and the motor can be shut down instantaneously 
by a magnetic disc brake. 

The end plug has a 3/4 inch threaded shaftway machined 
into it. Between the end plug and the shaftway is a .27 
inch hole to allow for the insertion of the pressure 
transducer. The stainless steel shaft is hollow to allow 
for the passage of the thermocouple and / the transducer leads 
to the slip ring unit. The shaft has 3/4 inch of its 
outside diameter threaded for screwing into the shaftway, 
and is threaded for 3/8 inch along its inside diameter on 
both ends. The condenser end mates with the transducer arm 
described earlier and the opposite end is joined to a 
coupling connected to the slip ring unit. 

There are two main bearings supporting the heat pipe. 

The large, 5.96 inch diameter, floating, single bearing 
seats on the evaporator flange which is supported by an 
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aluminum pillow block. Copper cooling lines through which 
tap water flows have been placed in the pillow block to 
help cool the bearing. There is an oil input line at the 
top of the pillow block which allows for lubrication of the 
bearing. The small double row, angular contact bearing on 
the condenser end is mounted on a seat machined on the shaft 
and is self-lubricated. 

b. Test Stand 

The test stand was designed and built so that 
the heat pipe could be rigidly supported with its longi- 
tudinal axis ranging in position from 0 to 90 degrees from 
the horizontal. Both the heat pipe and variable drive 
motor are mounted on steel bed-plates connected by a 4 inch 
iron pipe. The iron pipe is held in place by three 2-inch 
thick clamps which are supported two feet off the ground 
by a 1/4 inch steel plate. As the heat pipe is rotated in 
any position, the drive motor rotates with it and the entire 
assembly is then securely held in place. The motor has been 
vibration mounted in order to keep the vibration to a 
minimum. 

c. Spray Cooling Assembly 

Cooling of the condenser section is done by 
spraying tap water on the outside of the condenser during 
rotation. The condenser section is enclosed by a 9-inch 
long 13-inch diameter stainless steel cylinder, the con- 
denser box. The ends of the condenser box have Gar lock 
rubber seals separating them from the condenser section. 
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The condenser box is slid into place over the condenser 
section through its open ends and is then bolted in place 
to the support plate. An 8-inch x 6-inch window has been 
cut in the condenser box in order to be able to make 
adjustments or repairs without removing the condenser box. 
During operation, a steel plate is bolted over the window. 
The condenser box has eight spray nozzles mounted at the 
same axial position. The nozzles may be varied to give a 
range of droplet sizes from 300 to 600 microns. 

The coolant flows from copper feed lines to a mixing 
tube located at the top of the condenser box. A quartz 
thermometer records the temperature of the coolant in the 
tube. From the tube, the coolant flows through plastic 
tubing to each of the spray nozzles. The coolant drains 
from the bottom of the condenser box to another mixing tube 
where the temperature is recorded by the other quartz 
thermometer. From there the coolant flows to an external 
drain . 

d. Safety Shields 

During operation, the entire heat pipe assembly 
is surrounded by 1/8 inch thick stainless steel shielding 
bolted to the support plate to ensure safe operation. 

The evaporator shield, shaped like an inverted U, is 
bolted to a 1.0 inch thick plexi-glass window on its open 
end. This window serves as a safety shield and also allows 
the internal pipe operation to be observed. The other end 
mates with the condenser box. 
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B. EXPERIMENTAL PROCEDURES 

The experimental procedures consist of: (1) Filling 

Procedure and (2) Operational Procedure. 

1 . Filling Procedure 

In order to obtain film condensation in the conden- 
ser, the inside surface of the condenser must be prepared 
so that proper "wetting" will take place during operation. 
The following procedure is used: 

a. Scrub the condenser surface with 190 °F Alconox 
detergent solution. 

b. Flush with distilled water. 

c. Spray with alcohol. 

d. Rinse with distilled water. 

e. Spray with acetone. 

f. Rinse with distilled water. 

Having completed the above procedure, the heat pipe is 
ready to be filled. In this study, distilled water was used 
as the working fluid. 

In order to get a sufficient range of saturation tem- 
peratures and not exceed the 35 psia limit on the heat pipe, 
distilled water ice was placed in the evaporator and a low 
vacuum was drawn. The filling procedure was as follows. 

a. System was cooled by circulating cooling water 
through the condenser nozzles and by inserting distilled 
water ice in the evaporator. 

b. Six distilled water ice cubes (approximately 
150 ml) were subcooled by flashing them with liquid N 2 and 
placed into the evaporator with tongs. 
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c. Glass windows were put into place and the flange 
was tightened with a 30 in-lb. torque. 

d. The vacuum was drawn through the glass tube 
connected to a mechanical vacuum pump. 

e. While the vacuum was being pumped, the glass 
tube was sealed using a gas flame. 

2 . Operational Procedure 

Having completed the filling procedure, the heat 
pipe is ready to run at a pre-determined PPM. About 10 
drops of lubricating oil are added to the oil fill tube 
at the top of the large pillow block. The cooling water 
lines are turned on so that the rotameter is at approxi- 
mately 20% flow. 

The drive motor is turned on and the heat pipe is 
slowly brought up to the desired RPM. Once a liquid 
annulus has formed in the evaporator and the desired RPM 
has been reached, the DC motor generator is turned on and 
power is added to the evaporator. Having reached thermal 
equilibrium for a particular power setting, data was taken 
by the instrumentation described earlier. 

Shut-down is accomplished by first shutting off the 
current to the heater, and then decreasing the RPM by use 
of the remote control unit. 

C. EXPERIMENTAL RESULTS 

The heat pipe was run at 702 RPM and 1404 RPM. The run 
of 702 RPM was chosen because it was the lowest RPM at which 
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centrifugal forces dominated the gravitational effects. 

Power inputs to the heat pipe ranged from 1 kW to 9 kW. 

1 . Preliminary Results 

Using distilled water as the working fluid, the heat 
pipe was rotated at 702 RPM to obtain the heat flow through 
the evaporator wall for a particular power setting. The 
method used to compute / the heat flow through the 

evaporator wall, and the power supplied to the heater coils 
is shown in Appendix B. Table II shows a comparison between 
the electrical power into the heater coils and the heat flow 
through the evaporator wall for that particular power 
setting . 

The data taken from the four thermocouples in the 
evaporator wall is obviously incorrect since part of the 
data shows that the heat flow through the evaporator wall 
is higher than the electrical heat in for that particular 
power setting. In other cases it appears to be low. 

There are two explanations for the incorrect measure- 
ment of the heat flow through the evaporator wall. The 
first is the fact that the reading of thermocouple #2 
appeared erratic. When the heat pipe is not rotating, 
thermocouple #2 appears to be reading correctly, but when 
the heat pipe rotates it begins to read inadequately. It 
was concluded that the problem existed in the slip ring and 
may be due to insufficient mercury in the slip ring. The 
second explanation pertains to the reason why the heat flow 
through the evaporator wall is higher than the heat into 
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TABLE II 



POWER IN AND HEAT FLOW THROUGH EVAPORATOR WALL 



POWER IN 
(BTU/HR) 


HEAT FLOW IN 
EVAPORATOR WALL 
(BTU/HP) 


3,345 


4,059 


6,712 


10,315 


9,963 


12,848 


13,325 


7,610 


16,988 


9,126 


20,751 


24,041 


24,031 


13,437 


27,278 


31,986 



the heater coils for a particular power setting. It was 
concluded that the outside edge coils on each side of the 
heating surface were conducting axially. As a consequence, 
there wasn’t a uniform temperature distribution axially 
along the heating surface as assumed. Therefore, the actual 
heating surface length was less than the assumed length in 
the calculation of the heat flow and would account for the 
high values of the heat flow obtained in Table II. 

At this point, it was decided to obtain the heat flow 
out of the heat pipe by performing an energy balance on the 
condenser cooling section. The coolant flov; rate was 
determined by a flow rotameter, and its mean temperature was 
measured before and after spraying, by quartz thermometers 
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mounted in mixing tubes located above and below the conden- 
ser box. Having this information, along with the coolant's 
specific heat, the heat flow was determined by a simple 
energy balance. In order to make this method more accurate , 
the inside of the condenser box was insulated with silicone 
foam rubber and the plastic coolant lines were insulated 
with fiberglass. Appendix B contains a table showing the 
heat flow out of the condenser at different values of the 
heat pipe saturation temperature. 

2 . Visual Observations 

Using a strobe light visual observations were made 
through the view ports. At 702 RPM, the light provided by 
the strobe light was insufficient. To increase the light, 
the strobe light's frequency was increased to three times 
the RPM of the heat pipe. As a consequence, the observations 
made at 702 RPM were not clear. The observation of the 
liquid vapor interface in the evaporator at low saturation 
pressures indicated violent boiling with large vapor 
bubbles breaking through the liquid vapor interface at 
702 RPM. At higher saturation pressures, the boiling was 
more controlled with small bubbles constrained to the 
surface of the liquid. At the same saturation pressure, the 
boiling was less violent at 1404 RPM. 

The visual observation in the condenser indicated that 
there were ripples and random break-up of the condensate 
film, exposing parts of the condenser surface to the vapor. 
Furthermore, the condensate swirled and washed the condenser 
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wall at a higher RPM than the wall itself. It is doubtful 
that the above observed effects were due to condensate 
turbulence. According to Rohsenow and Choi 13], in the 
presence of interfacial shear, turbulent condensate flow 
occurs for a Reynolds number between 100 and 1800. For the 
low vapor velocities which occurred during the observed 
experimental runs , the calculated transition Reynolds number 
in the condensate was approximately 1000. However, the 
highest Re calculated for the condensate was only 291. It 
was concluded, therefore, that the swirling and random film 
break-up were not caused by condensate turbulence but were 
probably caused by vapor interaction with the condensate. 

Except for the 1 kW case, the Reynolds number calculated 
for the vapor flow for all of the power inputs was between 
2000 and 9000 (for 1 kW, Re = 972) . Since flow in a smooth 
tube is turbulent for Reynolds numbers between 2000 and 
50,000 [3], it was concluded that the vapor flow was tur- 
bulent except for the 1 kW case. The swirling of the 
condensate was probably caused by the higher velocity vapor 
striking the condensate. 
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III. DISCUSSION OF RFSULTS 



The graphs in Figures 9 and 10 show the experimental 
and theoretical heat transfer rates as a function of the 
saturation temperatures (corresponding to an internal vapor 
pressure) . 

In order to find the theoretical values of the heat 

transfer rate, the value h, the heat transfer coefficient 

was needed. The method used in determining h is shown in 

Appendix B. Because of experimental uncertainties, the 

values of h . = 1000 and h = 1800 were used since these 

mm max 

values were very close to the average of the minimum and 
maximum experimental values of h determined for the dif- 
ferent operating conditions. 

As observed in Figure 9, the theoretical values of Q m ^ n 

and Q are less than the experimental value of Q at 
m 3.x 

702 RPM. In the development of the analytical solution it 
was assumed that smooth film condensation occurred in the 
condenser. The visual observations described in the experi- 
mental results indicated that a smooth condensate film did 
not exist for the operating conditions used in these 
experiments but rather that the condensate film swirled and 
washed the condenser wall and appeared to break up in a 
random manner exposing part of the condenser surface to the 
vapor. The swirling effect and the resulting increased 
convection in the condensate increased the heat transfer 
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25.0 



Saturation pressure (psia) 
1.7 4.7 11.5 




Saturation temperature (°F) 

702 RPM 
Too = 67°F 
4 > = 1 ° 

Figure 9. Feat transfer rate vs saturation temperature 
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rate. It was also possible that the ripples and film 
break-up decreased the overall film thickness, reducing the 
thermal resistance, and thereby increasing the heat transfer 
rate. 

Figure 10 shows that the analytical and experimental 
results compare more favorably for the 1404 RPM run. There 
are two probable explanations. First of all, the higher 
RPM could have flattened the film causing it to approach a 
smoother film than in the 702 RPM run. If this were the 
case, the smooth film condensation assumptions would be 
better, and the analytical solution would be closer to the 
experimental results . Although the visual observations 
seemed to indicate that there were no non-condensibles 
present, the presence of non-condensibles in the condenser 
could be the second explanation. If non-condensibles were 
present, the experimental heat transfer curve would be 
shifted to the right since the heat transfer rate would 
require higher vapor pressures (and, therefore, higher 
saturation temperatures) . 

If the inlet coolant temperature were held constant, 
the heat transfer rate should increase with increased RPM 
for the same saturation temperature. For example, the 
analytical solution for at 1404 RPM at approximately 

the same used at 702 RPM resulted in higher heat transfer 
rates as shown by the dashed plot in Figure 10. Note that 
the experimental data at 1404 RPM does not appear to show 
an increase in heat transfer rates from the 702 RPM run. 
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Saturation pressure (psia) 
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Saturation temperature (°F) 

1404 RPM 
Too = 73.6 °F 
* = 1 ° 



Figure 10. Heat transfer rate vs saturation temperature. 



48 



This is due to the fact that the 1404 RPM run was conducted 
at a higher T ro (7° difference) than the 702 PPM run. This 
increase in T^ increased the saturation temperature by 
approximately the same amount such that the experimental 
heat transfer rate data points were shifted to the right. 

At the same saturation temperature, therefore, the heat 
transfer rate was lower for the higher T w even though it 
was at higher RPM. 
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IV. 



CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

1. The rotating non-capillary heat pipe is an effective 
heat transfer device. 

2. Smooth film condensation was not observed during the 
experimental runs. 

3. The approximate numerical solution conservatively 
predicts the heat transfer rate with a deviation of 18% at 
702 RPM. The deviation of the numerical solution decreases 
to 5% for 1404 RPM. 

4. The performance of the heat pipe is strongly 
dependent on the method used to cool the condenser. 

B . RECOMMENDAT ION S 

1. Run the heat pipe at higher RPM to determine 
whether increased RPM provides better agreement between the 
analytical and experimental results. 

2. Increase the instrumentation of the condenser to 
give more accurate data for the experimental heat transfer 
rate out of the condenser. 

3. Revise the analytical solution to include the ripples 
and film break-up in the condenser. 

4. Obtain experimental results for the effect of non- 
condensibles in the heat pipe. 

5. Use different methods of cooling the condenser wall 

by using different spray nozzles. 



6. Use a condenser section with a half-cone angle 
greater than 1° to obtain experimental data for the effect 
of larger cone angles on the heat transfer rate. 

7. Use working fluids with boiling points lower than 
water in order to find optimum heat transfer rates at lower 
saturation temperatures. 
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APPENDIX A 



CALIBRATION 

1. Calibration of Pressure Transducer . The pressure trans- 
ducer consists of a temperature-compensated, semiconductor 
strain gage Wheatstone bridge mounted on a 1/4-inch diameter 
diaphragm. It was calibrated against a Wallace-Tiernan 
standard pressure gauge (with a known accuracy of 0.03 ± 

0.01 psia) for pressures greater than atmospheric, and 
against a mercury manometer (with an uncertainty of ±.02 
psia) for pressures less than atmospheric. The transducer 
was calibrated at different operating temperatures by 
placing it within a controlled temperature oven with a 
quartz crystal thermometer (with an uncertainty better than 
±0.1°F). Reproducible data was obtained from 0-30 psia 

and from 65-222°F, and the transducer sensitivity was 
measured to be 3.54 ±0.06 ;mv/psia. Although the transducer 
Wheatstone bridge was temperature compensated, a non-linear 
zero shift with temperature was measured and incorporated 
into a series of calibration curves. 

2. Calibration of Thermocouples . Nine 30 gauge, Kapton 
insulated, copper-constantan thermocouples were carefully 
manufactured at the Naval Postgraduate School. They were 
calibrated up to 384°F using the same temperature controlled 
oven and quartz crystal thermometer mentioned above. In 
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addition , the following fixed point calibrations were made: 

a. Triple point of water (32.0°F) 

b. Boiling point of water (212. 0°F at 14.696 psia) 

c. . Melting point of tin (449. 4°F) 

d. Melting point of lead (621. 3°F) 

There were a total of nine calibration data points 
taken over the range of 32 -621.3°F. Using the NBS table 
of calibration data for copper-constantan thermocouples as 
a zero reference, a deviation curve of the error (mv) vs. 
the thermocouple reading (mv) was plotted. This curve (with 
an uncertainty of 0.1%) was then used as a correction 
factor for each data point taken during the experiment 
before entering the NBS table for the temperature. Recali- 
bration data was taken at a later time and the repeatability 
of all nine thermocouples was excellent. Using the boiling 
point of water and the melting point of tin, the deviation 
was within 0.005 mv of the previous data points. 

3. Calibration of Rotameter . The method used in the cali- 
bration of the flow rotameter consisted of using a Toledo 
scale to measure the mass of the coolant flowing into a 
bucket placed on the scale. As coolant flowed into the 
bucket, an arbitrary starting point on the scale was taken 
and a stop watch was used to measure the time it took for 
the weight to reach 20 pounds. This procedure was repeated 
several times for each per cent flow rate used on the 
rotameter. Data points were obtained for 14, 16, 18 and 
19.7 per cent flow rates. From this data a linear calibra- 
tion curve was obtained. 
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APPENDIX B 



DATA REDUCTION 



1. Calculation of Q EVAP 

The heat flow through the evaporator wall may be reduced 
to a one dimensional steady state problem if it is assumed 
that there is no axial heat conduction and steady state 
conditions exist. Under these conditions, the well known 
conduction equation in cylindrical coordinates may be 
written II] : 



or 



1 _d_ 
r dr 



dT 

k(T) r ^ 

Ql 



0 



k ( T) r 



dT 

dr 



B 



Constant 



(a) 



where : 



k (T) = k Q + m(T- T ) 



(b) 



A plot of k vs. T was made based on data taken from the 
Metals Handbook [5] for the thermal conductivity at different 
temperatures. From this plot the slope, m, and k Q were 
determined . 

Substituting equation (b) into equation (a) and inte- 
grating yields: 



F(T) 




B In r + c^ 
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Using a sub-routine for the method of least squares, 
the IBM 360 computer was used to find the slope, B, of 
F (T) vs. In r. From Fourier's law: 



■^EVAP = L 2tt 



[k(t) r || 



q evap 



L 2 tt B 



where 



L = the effective heating surface length in the 
evaporator 



EXAMPLE: T gAT = 153. 7°F, 702 RPM 



Q EVAP 



L 2 77 B 



L = .328 (ft. ) 

B = 4430. (BTU/hr - ft.) 

Q EVAP = 9,126 (BTU/hr) 



2 . Calculation of POWER IN 

P IN = Voltage * Current 



3 . Calculation of Experimental Q 



OUT 



°OUT m C p (T 1 " T 2 ) 



where 

m = mass flow rate (lbm/hr) 

Cp = specific heat (BTU/lbm °F) 

T^ = inlet coolant temperature (°F) 
T 2 = exit coolant temperature (°F) 
Q out = heat flow (BTU/hr) 



EXAMPLE: TSAT = 153.72°F, 702 RPM 



m = 797 ± 18 

C p = 1 *° 

(T 1 '" T 2 ) = 16 * 04 ± °* 31 
Q out = 12,784 ± 541 



4 . Calculation of h, outside condenser wall heat transfer 
coefficient 



Q OUT = hA c (T w " T j 



h Q OUT //A c (T w " 



where 

2 

A c = cooling surface area (ft ) 

T tt = outside wall temperature of condenser (°F) 
= coolant inlet temperature (°F) 
h = heat transfer coefficient (BTU/hr-f t 2 - 0 F) 



EXAMPLE: TSAT = 153.72°, 702 RPM 

Q qut = 12,784 ± 541 

A = .34 
c 

(T w - T^) = 28.44 ± 6.7 



h = 1405 ± 386 



5. Tables of Reduced Data 



TSAT 

(°F) 


m 

(LBM/HR) 


T - T 
1 2 

(°F) 


T 

w 

( 


- T 

oo 

°F) 


q out 

(BTU/HR) 


h AVG , 

(BTU/HR-FT -°F) 








702 


RPM 






93.4 


835± 23 


3.061.11 


6.611 1.8 


2,5551165 


12421411 


107.3 


832123 


5.921.25 


11.1 


1 2.6 


4,9251350 


14031423 


122.8 


804120 


9.231.31 


17.2 


1 2.7 


7,4211440 


13051280 


137.8 


823124 


12.5 1.54 


22.9 


1 4.6 


10,2961757 


13881375 


153.8 


797118 


16.0 +.31 


28.4 


1 6.7 


12,7841541 


14051386 


172.1 


764119 


20.5 1.25 


34.1 


1 7.2 


15,6471585 


14171350 


188.3 


764119 


23.6 1.31 


38.9 


110.3 


17,9921690 


15151465 


204.4 


863119 


22.6 ±.34 


37.9 


1 8.5 


19,5381730 


16491438 






] 


L404 : 


RPM 






98.4 


870123 


3.021.04 


1.0 


11. 1 


2,627+102 


12251408 


115.3 


850124 


5.961.18 


12.6 


15.2 


5,0661300 


14451664 


141.8 


825122 


12.2 1.21 


21.7 


. 0 


10 ,0 821433 


14.171 320 


174.3 


790+20 


18.9 1.53 


31.7 


16.5 


14.9391783 


1456 ± 371 


205.6 


786127 


25.4 1.41 


36.1 


13.9 


19 ,964 + 1011 


16431259 


224.3 


786127 


29.4 1.68 


45.6 


19.0 


23,13211347 


15631395 
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